Continuous Neel to Bloch Transition as Thickness Increases: Statics and Dynamics 
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We analyze the properties of Neel and Bloch domain walls as a function of film thickness h, for 
systems where, in addition to exchange, the dipole-dipole interaction must be included. The Neel to 
Bloch phase transition is found to be a second order transition at he, mediated by a single unstable 
mode that corresponds to oscillatory motion of the domain wall center. A uniform out-of-plane 
rf-field couples strongly to this critical mode only in the Neel phase. An analytical Landau theory 
shows that the critical mode frequency uo ~ ^Jhc ~ h just below the transition, as found numerically. 




PACS numbers: 75.60.Ch, 75.70.Ak 

Keywords: phase transition, domain wall, Neel, Bloch 

Introduction. Domain walls in ferromagnetic thin 
films, including exchange, uniaxial anisotropy, and the 
dipole-dipole interaction, have been extensively studied 
for the past 100 years. [1] For a monodomain of small 
film thickness h, the dipole-dipole interaction causes the 
magnetization to lie entirely in the plane. Such domains 
in a thin film are separated by either Bloch or Neel do- 
main wall. For the Bloch domain wall the transition be- 
tween domains occurs with the magnetization developing 
an out-of-plane component. [51 [3] The system develops 
surface poles, and the associated total dipole-dipole en- 
ergy is proportional to h. For the Neel domain wall the 
magnetization lies entirely in the plane of the film.jH |S] 
The total dipole-dipole energy for a Neel wall comes from 
the self-interaction of a volume pole density, and thus is 
proportional to . It is thus clear that for small thick- 
ness h the Neel wall is energetically preferable, whereas 
at larger h the Bloch wall wins. This Neel-Bloch domain 
wall transition is the major focus of this letter. 

Recent experimental and micromagnetic studies show 
that a Bloch wall has a complex structure. [51 [7] Its 
interior magnetization has an out-of-plane (Bloch-like) 
component, whereas its near-surface magnetization is in- 
plane (Neel- like) — the so-called Neel caps (Figure 1). 

The Neel domain wall, on the other hand, has a nar- 
row central part which resembles the usual "exchange" 
domain wall and long logarithmic tails |8] [9]. These tails 
are the consequence of the dipole-dipole interaction. 

When the sample is sufficiently wide the so-called 
cross-tie domain wall can be an alternative equilibrium 
configuration. [TO] In this letter we completely neglect 
such configurations. We assume that at the thicknesses 
of interest the system possesses a single domain wall. For 
h less than the critical thickness he, the domain wall is 
of Neel type, whereas for h > he the magnetization in 
the domain wall itself tips out-of-plane, forming so-called 
symmetric and asymmetric Bloch walls, [llj For the thick- 
nesses considered here, Neel walls and asymmetric Bloch 
walls are the only low energy states. 



Figure 1: Neel and asymmetric Bloch domain walls, magne- 
tization distribution in the x-z plane. Sample (permalloy) 
is infinite along y, has thickness h along z and has width 
w = 200 nm along x. 

Despite its fundamental nature, it has been unclear 
whether the Neel to Bloch transition at he is continuous 
or discontinuous. [121 1131 HI] A recent and comprehen- 
sive magnetostatics study by Kakay and HumphrevjlSj 
indicates that the transition between Neel walls and 
asymmetric Bloch walls is first order, although the au- 
thors were cautious in identifying the specific thickness 
at which the transition occurs. 

The present Letter considers the nature of this transi- 
tion (continuous or discontinuous) by extending our work 
on the statics[9] and dynamics for Neel and Bloch walls. 
We numerically study f 16' the spectrum of normal modes 
as a function of the sample thickness h, and identify 
an unstable mode whose frequency vanishes at a criti- 
cal thickness he- We numerically show that the ground 
state of samples with h < he is a, Neel domain wall and 
the ground state of samples with h > heis an asymmetric 
Bloch wall. 

We also show that the transition can be described ana- 
lytically by a Landau theory of second order phase tran- 
sitions with a single order parameter — the amplitude 
of the unstable mode. We derive the Landau free energy 
for the transition up to the fourth order in the order pa- 
rameter. In agreement with the numerics, the frequency 
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of the critical mode uj ^ \/hc ~ h just below the critical 
thickness. 

Micromagnetics. Consider a ferromagnetic strip of 
thickness h along z, width 2w along x, and infinite along 
y (inset to Figure 1). We assume that the exchange 
length lex — \/ A/iTrM"^ satisfies l^x <C w. Our sam- 
ples are large enough for the domain wall center to com- 
pletely fit within the sample. The material parameters 
chosen are appropriate to permalloy: exchange constant 
A = 1.30 X 10~® erg/cm and saturation magnetization 
Ms = 795 emu/cm'^, which gives lex = 5.72 nm. We take 
no crystalline anisotropy. The magnetization for x = 
is taken to be parallel to x-axis and free boundary con- 
ditions are taken on the sides of the strip. 

The calculation of the spectrum for samples with dif- 
ferent thicknesses started from thin samples h (1-5 nm) 
and magnetization m° = yMg, x < 0, rfif = —yMg, 
X > 0. We equilibrated and then calculated the normal 
modes and their coupling to a uniform rf external mag- 
netic field. Using this equilibrium as a starting point in 
the relaxation algorithm, we then gradually increased h 
and calculated the new equilibrium configuration, normal 
modes, and rf couplings. We also started with large thick- 
nesses and gradually decreased the thickness, to check for 
bi-stable solutions. 
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imately corresponds to a z-dependent oscillation of the 
domain wall core along x. The oscillations are confined 
to the central part of the Neel domain wall, whereas the 
logarithmic tails remain unperturbed. The frequency uj 
depends on both the width 2w and on the height h. 

To a very good approximation, the (un-normalized) 
critical mode has in-plane oscillations of the form 
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Here 5 - the half- width of the domain wall center - was 
determined from a fit to the central region of the domain 
wall. The x-y phase factor (/)(a;, z) determines the sym- 
metry properties of the mode, but is not needed for the 
present purpose. 

The out-of plane oscillations rn'z\x, z) have the form 
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where the z phase factor (l)z{x, z) differs from 4>{x, z) . The 
amplitude of the out-of-plane oscillations is typically less 
than 25% of the amplitude of the in-plane oscillations, 
except for h near he- 
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Figure 3: Magnetostatic energy and its first derivative as a 
function of tliickness h {w — 100 nm). 



Figure 2: Frequency lj of the critical mode as a function of 
sample thickness h for different half widths w. Inset: criti- 
cal mode frequency near the critical thickness he, numerical 
results (dots) and square root fit (line) for w — 100 nm. 

The main two results of the micromagnetics calcula- 
tions are: (1) The ground magnetization state for small 
thicknesses is indeed the Neel domain wall. In this sys- 
tem the Neel wall has three distinct regions — the central 
region of width 26 and logarithmic tails. |9: (2) There is 
one mode of oscillations with frequency to which goes to 
at some critical thickness he (Fig 2). This mode approx- 



The inset to Figure 2 shows, for w = 100 nm, the thick- 
ness dependence of uj just below the critical thickness he 
(dots) as well as a fit (line) confirming the ^/he — h de- 
pendence. As h increases above he, and ground state is 
taken as Neel wall, the imaginary part of iv goes from 
negative (stable) to positive (unstable). This is the only 
unstable mode. 

Near he conventional methods (conjugate gradient al- 
gorithm, Landau-Lifshitz equation in the presence of 
large damping, etc.) failed to find the equilibrium state. 
This is because among all the modes in the Neel state 
only one is unstable; therefore any random perturbation 
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of the initial configuration will contain only a negligible 
projection on the unstable mode. Moreover, the unstable 
imaginary part of the eigenfrequency is of order 1000 Hz, 
which would require an enormous time interval for the 
instability to be observed by numerical integration. The 
only method we found to produce reliable and consistent 
results was to study the normal modes for the unstable 
initial configuration, find the unstable mode, and then 
add a component of an unstable eigenvector to the static 
solution. This led to a significant decrease of energy and 
gave a magnetization close enough to the true local equi- 
librium that more conventional methods then applied. 

For w = 100 nm the domain wall for h > he, where 
he ~ 39 nm, is indeed an asymmetric Bloch wall (Fig- 
ure 1). Figure 3 shows that to numerical accuracy both 
the total magnetostatic energy and its derivative remains 
continuous through the transition, although the dipole- 
dipole and exchange energies individually have discontin- 
uous first derivatives. The second derivative of the mag- 
netostatic energy is slightly discontinuous at the transi- 
tion, which implies a second order transition. The dif- 
ference between this and previous work |14j may be due 
to the failure of conventional numerical methods at the 
transition. If the calculations are performed with a false 
equilibrium (metastability), then as the film gets thicker 
more modes become unstable. At some thickness the nu- 
merical method finally finds the true equilibrium; at this 
point the transition would appear to be sudden. 

We also numerically studied the dependence of he on 
the exchange length l^x = ^/AJ2nM^ and the sample 
width w. By varying these parameters and later fitting 
the results to about 5% we find the best fit to be given 
by the form (Figure 4): 
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Figure 4: Critical thickness he as function of width w. dots 
give numerical results; curve gives fit with Q. 



Analytics. This section studies the symmetry prop- 
erties of the transition mode and its corresponding state. 
It then relates the width i5 of the Neel domain wall center 
and the thickness of the strip h at the transition. 

For very thin strips {w ^ h) below the critical thick- 
ness he the stable magnetization configuration M^{r) is a 
Neel domain wall (Fig. lb) with center at a; = and par- 
allel to the z-y plane (Fig. la). Because h is very small 
is negligible, and M° and My are nearly independent 

of z. (Due to translation invariance M does not depend 
on y.) Thus the most general Neel domain wall is 



MO(r) = Msh = Ms{u{x),v{x),0), 



(5) 



where u{x)'^ + v{x)'^ — 1, subject to u{ztw) — 0, v{±w) — 
=F1, with u{—x) = u{x) and v{—x) = ~v{x). For i5 <C w, 
near the center of the domain wall, u{x) w sech^x/S) and 
v{x) « ~tanh{x/S). 

Above the critical thickness he the stable configuration 
is the asymmetric Bloch wall (Fig.lc), whose magneti- 
zation M^if) lacks the x and z reflection symmetries, 
but preserves inversion symmetry {x,z) (— x, — z), 
Mx,z{-x^ -z) = M^,zix, z), M^{~x, -z) = -M^(x, z). 

The critical mode responsible for the Neel-to- 
asymmetric-Bloch-wall transition has magnetization 
jl{x,z) = M^(x,z) — M'^{x,z) and is symmetric under 
reflections of x and z, so 

f^x{ X , Z^ /i-j,(x, z), flx,y{x, z) flx.yi-^T Z^ 

liy,z{-x,z) = iiy^z{x,z)\ ^iz{x,-z) = fizix,z). (6) 

Since the amplitude of the critical mode is small, its 
magnetization is orthogonal to the static magnetization 
and can be written as 

fi = Ms{X{x, z)v{x), ~X{x, z)u{x)X{x, z)), (7) 

where A(— x, z) = X{x,z), X{x,—z) = — A(x, z), Ci^iz) — 
(^{x, z) and C(a;, —z) = C,{x, z). 

Expanding in small z subject to the above symmetry 
conditions, with z — 2z/h, we use the approximations 

A(x, z) = (Aiz -I- A3z3) /(x), C,{x, z) = (Co + (2^') g{x), 

where Aq, A2, Cii ^^^d Ca are constants to be determined. 

The micromagnetic calculations show that the critical 
mode is localized near the center of the strip and that 
it shifts the Neel domain wall core parallel to x. Thus 
l^x{x,z) - i^M^{x) and Hy{x,z) ~ i^M^{x). By com- 
parison with the flts p]|3| to the micromagnetic calcula- 
tions, a good set of approximations is 

f{x) = sech {x/5) , g{x) — sech {x/5) [2sech {x/5) — 1] . 

(9) 

The mode energy now takes the form 

^] = Wex[n, + Wdd[n, m]; (10) 
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where Wex and Wdd are the exchange and dipole-dipole 
contributions. Although the mode is localized in the cen- 
tral part of the domain wall, the exchange contribution 
to its energy is negligible, so[T7] 

Wm[n,p] = I dridf2{iJ.{ri)dr^){iJL{r2)dr. 



(11) 

We now expand the free energy functional using the 
trial solution (7]|9l. The energy of the mode is calculated 
to second order in the Aj, and by stability has no first- 
order terms. With the row vector Z — (Ai, A3, Coj C2) and 
its column (transpose) vector Z"^, we have 



(12) 



where Wmode is a 4 x 4 matrix whose elements depend 
on the system parameters. 

At the critical thickness he one of the eigenvalues of 
Wmode goes to zero, so its determinant goes to zero: 



DetW„ 



0. 



(13) 



For small exchange length l^x, as we have here, Wex ^ 
Wdd- The term Wdd is a function only of S/h. Solving 
(flfl yields ~ OMc (to about 10%). (Here S ^ 6^ is 



evaluated at the transition.) This is in good agreement 
with the numerics, where for different w and lex we find 
0.74 < hc/Sc < 0.78. The unnormalized critical vector is 



Ze = (0.19,0.70,0.66,0.17). 



(14) 



This yields the form of fl(x,z); its amplitude is deter- 
mined by a higher order expansion in the energy. 

The smallest eigenvalue Aq of the matrix Wmode goes to 
zero linearly when h approaches he- Because the rest of 
the modes have finite energy a standard result of Landau- 
Lifshitz equations is that the frequency of the mode cj^ ~ 
Ao ~ /ic — h, so that the frequency goes to zero at the 
critical thickness as a square root u ~ \/he — h. 

The expression (7) gives only the first order approx- 
imation to the transition mode fl{x,z). The second or- 
der correction xi^i ^) can be obtained from |M^(x, z) \ — 
|M°(x, z) + fl{x, z) + x{x, z)\ — Ms and fl{x, z) _L x{x, z). 
Using this correction, one can calculate the free energy 
near the transition: 



Wmode{'n)=V^Q.iM^sLh' 



he 



O.OOSry^M^i/i^, 
(15) 

where -q is the amplitude of the critical mode. 

Because there is no rf term, the relation between 
Wmode{h) and the mode frequency lo is: 



V2 
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dWmode 
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h=hc+0 



Summary and Discussion. Experimental study of 
the transition by measurement of the magnetization and 



even the critical mode would be desirable. Because the 
symmetry under ~z to +z of the critical mode changes 
at the transition, from in-phase to out-of-phase, there 
are drastic symmetry changes in the absorption. Thus, 
whereas for h < he (Neel wall) the mode strongly couples 
to a uniform out-of-plane rf field, for h > he (asymmetric 
Bloch wall) it becomes anti-symmetric with respect to the 
z axis and cannot couple to this same rf field. 

To summarize, we have found that the Neel to asym- 
metric Bloch wall transition in dipole-dipole coupled thin 
magnetic films of infinite length along one direction, with 
relatively large width 2w, is a smooth function of thick- 
ness h. The normal mode frequencies reveal that at a 
critical thickness he the frequency of the transition mode, 
where the nature of the domain wall changes, goes to zero 
as \/he — h. To high accuracy, both static and dynamic 
calculations indicate that this transition is continuous, as 
also supported by analytical studies. 
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